Hexadecapolar Kondo effect in URu 2 Si2? 
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We derive the coupling of a localized hexadecapolar mode to conduction electrons in tetragonal 
symmetry. The derivation can be easily adapted to arbitrary multipoles in arbitrary environment. 
We relate our model to the two-channel Kondo (2CK) model and show that for an /^configuration, a 
relevant crystal field splitting in addition to the 2CK interaction is intrinsic to tetragonal symmetry. 
We discuss possible realizations of a hexadecapolar Kondo effect in URu2Si2- Solving our model we 
find good agreement with susceptibility and specific heat measurements in Thi_ ;E U a ;Ru2Si2 (a; <C 1). 

PACS numbers: 71.10.Hf, 71.27.+a, 75.20. Hr 



Introduction. In an seminal paper [1], Cox addressed 
important differences between U- and Ce-based heavy 
electron materials in terms of the atomic structure of 
their /-shell. In Ce-based systems the most probable 
configuration has one /-electron, in contrast to many U- 
based materials having an / 2 -many body state as the 
most probable one. In crystal structures where the U- 
site has cubic symmetry, the / 2 -states can give rise to 
quadrupolar degrees of freedom which when coupled to 
conduction electrons lead to two-channel Kondo non- 
Fermi liquid behavior [1, 2]. Meanwhile various multi- 
polar orderings have been observed [3] as well as pro- 
posed as candidates for "hidden order" (HO) in materi- 
als with clear phase transitions but without an obvious 
order parameter. A prominent example in this area is 
provided by URu2Si2- For this material quadrupolar [4], 
octupolar [5], hexadecapolar [6] and triakontadipolar [7] 
order parameters have all been put forth. Recent exper- 
iments have definitely ruled out quadrupolar order [8], 
whereas another might implicitly hint at it [9]. On the 
other hand, the hypothesis of active Uranium hexadecap- 
olar degrees of freedom provides a natural explanation 
for numerous experiments [10]. In this Letter, we gener- 
alize the work of ref. [1] and present a simple construc- 
tion of low-energy Hamiltonians that describe the cou- 
pling between multipoles and conduction electrons in the 
tetragonal crystal field of URu2Si2- We show that the U 
hexadecapolar degrees of freedom, couple symmetrically 
to multiple channels of conduction electrons. We solve 
the resulting model using the numerical renormalization 
group, and thereby successfully describe the properties of 
Th i _ 2; U a; Ru2Si2 (x < 1) [12, 13]. Thus hexadecapolar 
fluctuations also serve as an explanation for the anoma- 
lies observed Th 1 _ a ,U :! ;Ru2Si2 (x <C 1). 

Construction. To construct a tractable model, valid 
at very low-energies, we take into account only the 
two lowest-lying /-configurations with double occupancy. 
This is motivated by recent LDA+DMFT calculations 
which indicate that, while U has mixed valence in this 
material with /-electron occupancy between 2 and 3, 
two crystal field singlets with double occupancy and 
different symmetries have the highest probability [6]. 



The Uranium degrees of freedom are then described 
by a J = 4 multiplet, split by the tetragonal crystal 
structure. The ground state and the nearest excited 
level are, respectively, time-reversal and parity even, 
A 2g and Ai g basis states of the group & x D4/, 
with & — {1, T} the group of time-reversal, 1 
the identity, T the time-reversal operator; and D^h 
the tetragonal point group including parity. Viz., the 
lowest-lying singlets are: \A-2.g) = r^g (|4) — |— 4)), 
and \A lg ) = £5^ (|4> + | - 4)) + sin0|O), given in 

terms of the eigenvectors, | J z ) , of the operator J z in 
the J = 4 multiplet with the quantization axis cho- 
sen parallel to the c-axis of the crystal. To keep the 
equations short, we follow refs. [1, 2] and assume, the 
/-shell of the U atom hybridizes predominantly with 
I = 3, J = I conduction electrons. The conduction 
electrons at the local site, ipjjj can be classified into 
the four double-valued or spinor irreducible representa- 
tions (irreps), Tq p , V j p of the tetragonal double point 
group, Dih, with p — g/u for parity even/odd ir- 
reps (i.e. for I even/odd). Under time-reversal symmetry 
T ipjjjT- 1 = (-) l ~ J + J ' ^\j,_ Jz) ■ We set up a basis 
so that the a = ± components of the Kramers dou- 
blets, ^ (n) (where n enumerates doublets of the same 

type within the same J multiplet, and j = 6,7) comply 
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ing the same for annihilation operators. For the local 
conduction electron basis, we choose the following two 
independent r -j u Kramers doublets for creation opera- 
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Kondo Hamiltonians are made up of spin-flip and 
diagonal processes: Hk = H± + H z ■ When con- 
structing these two parts connecting the two singlets, 
the only relevant, non-trivial tensor products of irreps 
are r 6u <g> T 6u = T 7u ® T 7u = A lg © A 2g E g and 
A 2g ®A 2g = Aig [15]. Inserting the appropriate tetrag- 
onal Clebsch-Gordan coefficients, symmetry thus binds 
the most general form for the spin-flip and diagonal parts 
to be [15] 
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The couplings are real and must satisfy J 1 2 1 = J 2 1 % 
to ensure hermiticity, but otherwise arbitrary. We omit- 
ted processes including r 6 electrons as they decouple 
from the impurity. The hexadecapolar, i.e. "spin-flip" 
fluctuations are thus coupled to four species of con- 
duction electrons, namely to the two independent Y 7u 
Kramers doublets. T-L± has the structure of the two- 
channel Kondo (2CK) model where the role of spin in- 
dex is played by the index that distinguishes the two 
different I^'s, and the channels are distinguished by 
the Kramers indices. To make this correspondence 

more explicit, we introduce the operators • 
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transformation: 
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and perform the unitary 



-> \A'- 
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which al- 



lows us to rewrite T-L± in the standard notation 
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v\v<rZvV q uS +h.c.+OS x 



with q £ {a, b} , fi, v e {f, 1} . Here and in the fol- 
lowing, repeated channel (g) and spin (/i, u) indices are 
to be summed over; O = contains only conduction 
electrons [25]; er + = a x + ia v is composed of Pauli 
matrices; and S+ = S x + iS y = \A 2g )(A' lg \; S z = 

(\A 2g ){A 2g \-\A' lg ){A' lg \)/2;S- = S+ f ; 11 = 
\A 2g }(A 2g \ + \A' lg ){A' lg \. 

Discussion. Channel symmetry is a consequence of 
time-reversal symmetry. The operator, O S x is irrel- 
evant around the 2CK fixed point (and marginal in the 
free fermion scaling regime), as shown either by NRG cal- 
culations or using conformal field theory results [14, 16]. 
It does not destroy the 2CK state, since it neither breaks 
channel symmetry, nor lifts the spin degeneracy. Thus we 
must have 2 ^ J 2 1 in order for overscreening to oc- 
cur. This asymmetry comes up naturally e.g. if we start 



off with a spherical symmetric Anderson Hamiltonian, 
perform the Schrieffer- Wolff transformation to arrive at 
a Kondo- type of interaction and then project to the crys- 
tal field states, |-A 2s ), \Ai g ) at strong spin-orbit (i.e. jj) 
coupling, as described in refs. [1, 2, 14]. 

The diagonal part, H z cannot lead to non-Fermi liq- 
uid behavior by itself, but it can quite possibly destroy 
it. Channel symmetry is preserved by time-reversal sym- 
metry. However, the level degeneracies are lifted by the 
crystal field, both between the \4A 2 /i g ) states and also 

(n) 

in each screening channel between T 7u electrons with 
different n's. The dangerous terms are 
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Both types of crystal field splittings are relevant around 
the 2CK fixed point with scaling dimension i [14, 
16], and present in % z with the amplitudes Aj TOp = 
(j z 221 + J 111 -J z 222 -J z 112 )/2, and A cond = 
(J 221 - J z lll + J 222 ~ J z 112 )/4. In fact, they are 
the only possibilities for relevant perturbations, if chan- 
nel symmetry is intact [16]. 

Thus for this model to exhibit 2CK scaling in some 
temperature range, Aj mp and A conc i must fall below the 
Kondo scale, Tx- This necessarily requires fine-tuning, 
and the basic assumption of the A 2g — A\ g scenario — 
and, as we show below, of any other doublet-ground state 
scenarios — is that this accidental degeneracy is responsi- 
ble for the unique behavior of URu2Si2 among the large 
number of U-based heavy fermions. LDA+DMFT calcu- 
lations for URri2Si2 are indeed consistent with this acci- 
dental degeneracy on the scale of Tk [6] . 

A local, z-directed magnetic field results in the follow- 
ing leading additions to the Hamiltonian 



/i B i (S + - S ) 



H. magi ^ 
Hmagc OC /J.B 

J,e{±f,±f} 



for the two-singlet part and for the local conduction 
electrons, respectively. These terms have similar effect 
as H r z f /c . Namely, the impurity part, W z f + H magi , 
amounts to an effective magnetic field (or crystal field 
splitting) pointing into other than the z-direction. The 
same holds true for the conduction electrons with the ef- 
fective magnetic field/crystal field splitting being differ- 
ent in the two channels. Thus while Hk is not identical 
to the 2CK Hamiltonian, it flows to the same fixed point 
when the relevant perturbations, which split apart the 
two different I\ irreps, or the two local singlets, vanish; 
and the application of magnetic field thus breaks both 
the channel and the spin symmetry of the 2CK model. 

Comparison with Other Scenarios. For the order of 
crystal field levels in URu2Si2, other scenarios have also 
been put forward in the literature. In turns out that the 
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FIG. 1: (a) : Symbols: Molar susceptibility, x™f, of the 5/ 
electrons in Thi-xUxRusSia at x = 0.03 as the function of 
temperature in magnetic fields between H = and 5 T (H \\ 
c). (a) — (b) : Curves: The local, dynamic susceptibility of the 
two-channel Kondo model in the presence of magnetic field, 
computed with density matrix-NRG at T = . (c) : The 
crossover scale, Tb (defined as the intersect of the low-w and 
high-oj asymptotes of x( w )> see plot (&)) shows quadratic and 
linear //-dependence in the 2CK and local moment scaling 
regimes, respectively in good agreement with experiment; Tk , 
the crossover scale between the local moment and 2CK scaling 
regimes, is ~ 1.3 K. Notice that the experimental data for 
H = T can be equally well fitted with every NRG curve 
where H < 0.2 re- 



structure of the effective low-energy Hamiltonian for any 
two quasi-degenerate states is rather similar to that of 
the two-singlet case considered above. This applies also 
to a proposed E g (or T 5 ) doublet ground state [12, 17- 
19], formed by \E g x) 1 \E g y). Fluctuations within this 
doublet can couple to two products of irreps : to F7 ® T 7 
and also to r7(g)r 6 . Importantly, the two different irreps 
play again the role of spin, while their Kramers indices, 
connected by time-reversal, play again the role of channel 
index in the 2CK language. If the two lower-lying irreps 
were degenerate, the spin symmetry of the 2CK model 
would be unbroken and 2CK scaling would occur. Since 
the degeneracy is approximate, the system will eventually 
flow to a Fermi liquid fixed point. 

In the HO phase of URu 2 Si 2 , the ref. [6] 
proposed order parameter, (|^-2g) (A i g |) , is non- 
vanishing due to its real part which, in highest or- 



der of the multipole expansion, contains the ex- 
pectation value of the hexadecapolar A 2g tensor: 

'{Jl - J*)(jjy + JyX) + (jJy + JyJ X )(M ~ ^)] [6]. 

However, the same reasoning can be repeated for 
(\E g x)(E g y\), whose real part also contains this hexade- 
capolar ordering. These points are substantiated by the 
construction of the Hamiltonian for the E g ground state 
in the supplementary material section. 

Comparison with Experiment. It has long been rec- 
ognized that Xa the (magnetic) dipole susceptibility of 
Thi_ 2; U 2 ;Ru2Si2 along the c axis shows logT behavior 
at low-T (see refs. [12, 19] and Fig. 1) in accord with 
the 2CK descriptions corresponding to both scenarios. 
However, susceptibility and resistivity measurements find 
that the magnetic field (H) induced crossover scale to a 
Fermi liquid depends on H linearly, i.e. Tjj <x H v with 
77 = 1, which does not agree with the 77 = 2 behavior 
corresponding to the 2CK scaling regime [13]. To make 
contact with these experiments, we solved the model, Eq. 
(1) by NRG, and confirmed that it indeed flows to the 
2CK fixed point where OS x is irrelevant. Then we added 
a magnetic field, mimicked only by Eq. (3), to the 2CK 
model, and solved this model using an upgraded version 
of our density matrix-NRG code detailed in ref. [20] . The 
values of the magnetic field and the Kondo coupling were 
adjusted to fit the experimental data of refs. [12, 13]. 

Invoking uj/T scaling, we fitted the T-dependence of Xc 
by the dynamic susceptibility of the 2CK model in mag- 
netic field, as we trust our dynamic correlation functions 
(produced by the density matrix algorithm at T — 0) bet- 
ter than the thermodynamic quantities. Fig. 1 shows con- 
vincing agreement between theory and experiment apart 
from the small discrepancy for T > 30 K, i.e. for large 
energies where the resolution of NRG is limited. We ob- 
tained Tr « 1.3 K from the fit (see the caption of Fig. 1 
for further details on Tk)- This finding places the mea- 
surements in magnetic fields around the crossover region 
between the local moment and 2CK scaling regimes. In 
both regimes, scale invariance entails the hyperscaling re- 
lation, rj + v = 2 with v the critical exponent defined 
by x « H~ v . Thus for H = 0, the observed v = 
gives rj = 2, meaning that for T between 0.1 and 10 K, 
the system is in the 2CK scaling regime. In contrast, for 
H = 1 to 5 T, the experiments measure 77 = 1 resulting 
v = 1. Thus, we conclude, these magnetic fields in addi- 
tion to the ubiquitous, relevant crystal field splitting, are 
(slightly) larger than Tk and the system flows directly 
from the local moment regime to a one-channel Kondo 
fixed point without traversing the 2CK scaling regime. 
In Fig. 1(a), the ratios of magnetic fields to Tr-, fitting 
the susceptibility, further illustrate this point. 

By taking a closer look at the specific heat coeffi- 
cient in Fig. 2(a), we can reinforce these statements, 
and get another estimate for Tk- Two regimes for the 
given magnetic field values are clearly visible: The curves 
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FIG. 2: (a) : The 5/ electronic specific heat coefficient in 
units of 7? (C p , 5f /TR) of Thi_ a; U I Ru2Si2 at x = 0.07, in 
magnetic fields between and 12 T, reproduced from ref. 
[19]. (6) : 2CK specific heat coefficient from NRG for dif- 
ferent magnetic field values. The two sets of curves display 
the same trend. However there is a factor of ~ 4 difference 
in the magnitudes (depending on the precise value of Tk). 
We ascribe this difference to experimental inaccuracy, as we 
observe a factor of ~ 2 difference between the magnitudes of 
C Pl sf/T for the same material published in refs. [12] and [19], 
respectively. 



for H = 0.5 and 1 T slightly overshoot the curve at 
H = T and low-T, in contrast to the curves for H > 2 T 
which exhibit a bump at around T w H . The rise of 
7 = Cp^f/T for low- fields at low-T is reminiscent of 
the 2CK scaling regime, except that the measured In- 
dependence of 7 at low-T for H = is not quite logarith- 
mic. These observations can be explained by the pres- 
ence of an effective crystal field splitting, < A < Tk, 
already at zero magnetic field, and placing Tk between 
1 and 2 T. From the susceptibility fit we estimate : A is 
anywhere below about 0.2 Tk (c.f. Fig. 1(6)). These as- 
sertions are further confirmed by our NRG calculations 
for the specific heat coefficient (see Fig. 2(6) or Fig. 4 in 
ref. [21]). For H > T K , the bumps at around T w H 
correspond to the Schottky anomaly due to the Zeeman 
splitting between the two local states. 

In the 2CK model, if A < Tk , there is a non-Fermi 
liquid (NFL) region over 2 log (^) decades, since the 
splitting induced crossover scale to a Fermi liquid de- 
pends on A quadratically. Susceptibility measurements 
find an NFL region over at least one decade, putting an 
upper bound on the ratio, A/Tr- < 0.6 and giving the 
conservative estimate: Tk < 5 T and A < 3 T. 

Conclusion. Motivated by recent findings on the elec- 
tronic structure of URu2Si2 [6], we derived the Kondo 
coupling of localized hexadecapolar fluctuations to con- 
duction electrons in tetragonal crystal field. The deriva- 
tion can easily be adapted to arbitrary situations. The 
coupling has the form of the 2CK model plus relevant, 
spin symmetry breaking perturbations present even in 
magnetic zero field. The local degrees of freedom are 
symmetrically coupled to two different irreps of conduc- 
tion electrons. Solving the model with NRG we showed 
that the hypothesis introduced in ref. [6] to describe the 



hidden order in URu2Si2 can consistently account for the 
behavior of URu2Si2 in the dilute limit, Thi_ a: U a ;Ru2Si2 
(x <C 1). Nonetheless this behavior, while confirms that 
the material is strongly correlated, does not discriminate 
between different competing scenarios for the ground 
state-excited state sequence in this material. Since the 
two irreps of conduction electrons are not connected by 
symmetry, we expect a Fermi liquid to emerge at suffi- 
ciently low-energies. Hence the intermediate non-Fermi 
liquid regime, observed in URu 2 Si2, is a result of an acci- 
dental degeneracy. This accidental degeneracy is respon- 
sible for the unique properties of this compound in the 
dilute and dense limits, among the hundreds of known U- 
based heavy fermion materials. We found that the scale 
of the crystal field splitting and the Kondo temperature 
is smaller in Thi_ 2 U 2; Ru2Si2 (x <C 1) than in URu2Si2- 
The splitting between the two T7 irreps should be sensi- 
tive to the conduction electron filling, and we expect it to 
be larger in Lai_ a; U 2 ;Ru2Si2 [22] where clear Fermi liq- 
uid behavior is observed. In this context it is also worth 
pointing out that in Thi_ :r U a ;Ru2Si2 (x <C 1), the resis- 
tivity follows an approximate log T behavior with a nega- 
tive coefficient suggesting that a crossover to a Fermi liq- 
uid behavior takes place at sufficiently low temperatures. 
The study of the resistivity, however, will likely require a 
more realistic model for the diluted URu 2 Si2, including 
all bands present in the solid; also a more sophisticated 
approach to calculating the resistivity in non- Fermi liquid 
quantum impurity models than the ones present in the 
literature [23]; and corrections to the very dilute limit as 
was done recently in the context of the NMR experiments 
[24]. 
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SUPPLEMENTARY MATERIAL FOR "HEXADECAPOLAR KONDO EFFECT IN URU 2 SI 2 ?" 
Kondo Hamiltonian for a Tetragonal E g (or Ts) Doublet Ground State 

The structure of the Hamiltonian for two singlets, \A2 g ), \A\ g ) and for any other two quasi-degenerate states has 
many common features. The proposed E g (or T 5 ) doublet for the U ion ground state in Tli 1 _ 2 ,U 2 ;Ru2Si2 (x -C 1) 
is claimed to provide full symmetry protection for the 2CK state [12, 17-19]. We show that crystal field fine-tuning 
is needed in this case as well, to hit the 2CK scaling regime. On the vector space spanned by the E g ground state 
doublet 

3>-|3»], 
3> + |3»], 

we have four independent operators defined e.g. as il e = \E g x){E g x\ + \E g y){E g y\ 1 a v E = 
i(\E g x)(E g y\ - \E g y)(E g x\) , a E = \E g x)(E g x\ - \E g y)(E g y\ and of. = \E g x)(E g y\ + \E g y)(E g x\. They 
transform as A\ g , A2 gi B\ g and B 2g tensors of D^h , respectively [15]. In the Kondo limit, these operators can 
appear in the Hamiltonian only in certain combinations with the conduction electrons [14]. Namely, the B\ g and 
B 2g tensors can only couple to products of T 6 and T 7 electrons following the rule: T 6u ® Ti u — B lg © B 2g ® E g , 
i.e. 

m— 1 ^ ' n , m — 1 ^ 

+ Jxa 2 (*^> + Sr«_ + *If>- 5 if> + ) °i + h ^ 

Tn=l ' n,m— 1 ^ 

with arbitrary real couplings. Again, each allowed term describes a 2CK screening process, and channel symmetry 
is guaranteed by time-reversal symmetry One apparent difference from the two-singlet case is that Tq u conduction 
electrons combined with I^'s can also participate in real 2CK screening processes. Besides there remains also the 
possibility of the two I^'s forming the two screening channels. Once we single out the two lowest-lying channels, the 
two cases become very similar in that the diagonal processes for the E g ground state also include a relevant splitting 
between the different crystal field channels of conduction electrons. Thus an E g doublet ground state does not enjoy 
level degeneracy protection, and in this respect, it is not distinguished from other, two-singlet ground state scenarios 
in tetragonal symmetry. Full protection of the 2CK state can be achieved in a setting were cubic symmetry merges 
the Tq and Tj representation into a Tg quartet as in Cox's original proposal [1]. 
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